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Abstract 

Time reversal of quantum dynamics can be achieved by a global change of the Ifamiltonian sign (a hasty Loschmidt 
daemon), as in the Loschmidt Echo experiments in NMR, or by a local but persistent procedure (a stubborn 
daemon) as in the Time Reversal Mirror (TRM) used in ultrasound acoustics. While the first is limited by chaos 
and disorder, the last procedure seems to benefit from it. As a first step to quantify such stability we develop a 
procedure, the Perfect Inverse Filter (PIF), that accounts for memory effects, and we apply it to a system of coupled 
oscillators. In order to ensure a many-body dynamics numerically intrinsically reversible, we develop an algorithm, 
the pair partitioning, based on the Trotter strategy used for quantum dynamics. We analyze situations where the 
PIF gives substantial improvements over the TRM. 



In the last years, the group of M. Fink in Paris devel- 
oped an experimental technique called Time Reversal 
Mirror (TRM) that allows the time reversal of acoustic 
excitations [1]. An ultrasonic pulse, produced inside a 
control region (also called cavity) where it suffers mul- 
tiple scattering processes, is detected by several mi- 
crophones as it escapes through the boundaries. These 
transducers can also act as loudspeakers and the reg- 
istered signal is played back in the time reversed se- 
quence. Thus, the signal focalizes in the source point 
forming a Loschmidt Echo [2]. According to the exist- 
ing theory, an exact control of the wave function in the 
cavity would require the control of the wave function 
and the normal derivative at the boundaries. However, 
the reversal is quite good even when these conditions 
are not fulfilled by the experiment: the detectors might 
not enclose the cavity or the recording time period 
could be reduced to a fraction. Another surprising fea- 
ture of this time reversion procedure is that it shows a 
much better stability in inhomogeneus media as com- 
pared to ordered ones. This leads to numerous appli- 



cations in medical physics [3] and communications [4] . 
A first step to asses the errors is to develop a proce- 
dure that could achieve perfect reversal. This task was 
developed for the domain of quantum waves and was 
named Perfect Inverse Fiher (PIF) [5]. The PIF pro- 
cedure assures the exact reversion by injecting a wave 
function that compensates precisely the feedback ef- 
fects through a frequency dependent renormalization 
that involves the exact Green's function at the injec- 
tion sites. Here, we use a simple microscopic model 
that presents wave behavior and describes energy dis- 
sipation to show that the PIF procedure also applies in 
the classical domain. The model, represented in Fig. 1, 
is a variation of that of Rubin [6] : a surface oscillator 
with mass mo and natural frequency (jjq is coupled to 
a semi-infinite harmonic chain of bulk oscillators with 
mass m. We are interested in the time reversion of the 
initial condition where all the oscillators are in their 
equilibrium positions except for the surface one. The 
energy stored in the surface oscillator is expected to 
decay due to the effective friction produced by the "en- 
vironment" of light masses. The Hamiltonian is: 
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Fig. 1. Scheme of the model: a simple pendulum (surface os- 
cillator) is coupled to the bulk masses. 
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where uji is the natural frequency of the ith oscillator 
and Ui describes the displacement from the equilib- 
rium. Springs with elastic constant accounts for 
the coupling between first neighbors. For the proposed 
model, all springs and masses are equal and = 0, 
except for the heavier surface mass that suffers an ad- 
ditional haxmonic restitutive force (wo is finite). The 
exchange frequency Wx = \fKjm and the ratio a = 
m/mo < 1 have been chosen to set the system in the 
extended-oscillatory dynamical regime [7]. The equa- 
tions of motion in the frequency domain can be written 
in the matrix form 



^\lu)u{uj) = (cj^I - M) u{u}) = 0, 



(2) 



where D(u;) = (w^I — M) ^ is the resolvent associated 
to the dynamical matrix in the site basis 
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The resolvent provides the solutions to Eq. 2 with im- 
pulsive forces fi{t) = miAui{0)d{t). This is, 

^jit) = ^e-'-*Dj,i{uj)Aui{0) (4) 



2n 



(5) 



Notice that Dj^t) relates the jth displacement ampli- 
tude due to an initial condition of velocity in the ith 
oscillator. In general, the solution of the Eq. 2 in pres- 
ence of forces Fi{t) results 
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and can be rewritten as 



where Fi{u}) = fi{u;) + gi{ijj) is the Fourier transform 
of the force applied at mass i that is a sum of two com- 
ponents: an impulsive force and a shifting force. This 
last is able to produce an "instantaneous" shift A«j (0) 
in the position without changing its momentum. This 
would require that a first "impulsive kick" be followed 
by a compensating one: 

gi{t) = lim miAMi(0)i [6 {t + ^t) - S {t - fr)] , 

which in frequency domain means: 



gi{iL>) = — iwmiAwi(O). 



(7) 



Thus, in the time domain 



(8) 
(9) 



which serves as a definition for the position-position 
response function, also called the Green's function: 

Gj,i(w) = —luiDj^uj) = ioJXj,ii'^)^i- (10) 

On the other hand we will use linearity to write the 

observed displacement iu terms of Sui{t), the forced 
position shift accumulated in the unit time: 



Uj{t) 



Gj,i{t-t')8ui{t')At'. (11) 



We seek the injection function 5ui{t) that produces the 
exact reversion of the original wave within the control 
region, i.e. uf^{t) = u*{2tR - t) for tR<t< 2tR. Ac- 
cording to Ref. [5] , the perfect time reversal is possible 
if the dynamics starts and ends up without any exci- 
tation inside the cavity. Even when our system starts 
with an "excited" cavity, the lack of momentum at each 
mass ensures that forward and backwards evolutions 
arc identical. Once the decay signal is registered at the 
transducer for positive time, the earlier time values 
(corresponding to a fictitious injection) are also known 
and we build the function to be inverted at time tn 

, , [ <{tB^-t), -tR<t<0 

[Usit-t„), 0<t<tR 

and the injection function can be obtained in the fre- 
quency domain by 



(13) 



This equation defines the Perfect Inverse Filter (PIF) 
for a classical wave, where the injection prescribed 
by the TRM procedure appears now corrected by the 
Green's function Gs,s{t). The time evolution of the dis- 
placement amplitude is shown in Fig. 2 for several situ- 
ations. The recording time Ir is longer than the decay 
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Fig. 2. Displacement amplitudes. Up: initial condition at the 
surface oscillator. Down - left: registered signal in Xs- Down 
- right: response function in Xg. The chosen parameters arc: 
a = 0.25, 0)0 = l.SoJx, tn = lOOOoj^^ and = 10. 

time in the whole cavity, i.e. all the masses in the con- 
trol region have enough time to recover their equilib- 
rium positions. The injection functions of both proce- 
dures differs near the bandedges {uj — ±2ljx) indicating 
that TRM corrections are importants in cases where 
the whole spectrum is involved (typically in broad- 
band experiments). When all the masses in the cav- 
ity are in their equilibrium position, the injection at 
the source point Xs produces oscillations that propa- 
gate both sides of Xa . Only dynamics inside the cavity 
is fully reversed. Fig. 3 shows the local energy of the 
surface oscillator 

its decay and recovering presents the three temporal 
domains [8]: It begins with a quadratic dependence, 
continues with an exponential decay associated with 
a Self Consistent Fermi Golden Rule, and becomes an 
inverse power law when the energy return is compa- 
rable to the residual energy in the surface mode. The 
time reversed signal reproduces all the three regimes, 
even when they involve very different signal intensities. 
We attempt to asses the quality of the reversal by the 
Loschmidt Echo [2] of a wave function ip{t) normalized 
in the cavity: 

M{t)^\{ip,,,{t)\ip{2tii-t))f tR<t<2tR. (14) 

First, the inner product uses a metric tensor given 
by the Hamiltonian, i.e. we refer the recovered energy 
to that originally contained within the control region. 
Considering the case where tn = 1000wr\ the PIF 
yields Mpif = 0.999 while the TRM gives A'/trm = 

0. 982. This result is not representative of the evidenced 
by Fig. 4. Alternatively, we use the Euclidean met- 
ric tensor, obtaining Mpip = 1 and A/trm = 0.765, 

1. e. in PIF all the initial condition have been recov- 



Fig. 3. Decay and recovering by PIF procedure of the lo- 
cal energy in logarithmic scale. The insets show in detail the 
quadratic and exponential regimes. 
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Fig. 4. Recovering of the initial condition for the two proce- 
dures. 



ered whereas in TRM there is a spreading of displace- 
ments and velocities along the cavity. The reversed lo- 
cal energy in both procedures is compared in Fig. 5. 
The PIF procedure can not be distinguished from the 
ideal reversal. The TRM has two failures: the ampli- 
tude of the local energy at t = 2tii is always less than 
the initial case and the temporal regimes are delayed. 
In summary, the time reversion of the dynamics of a 
pendulum coupled to an harmonic chain has been done 
by means of two types of "stubborn daemons": 1) The 
TRM, which neglects memory effects (meaning that an 
instantaneous system response is assummed). 2) The 
PIF, that accounts for memory and feedback. This last 
is much better if the initial condition is build up over 
the whole range of frequencies. In this case, the correc- 
tion of the Green's function becomes non trivial and 
ensures a better reversion quality. 
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Fig. 5. Focalization of tlie local energy in logaritmic scale. Solid 
line represents the ideal time reversed decay, dashed line is PIF 
and TRM is shown in dotted line. Note that PIF is superposed 
with the ideal case. 

Appendix A. Classical dynamics through pair 
partitioning 



While the ID dynamics in the model considered 
can be obtained analytically by the continued fraction 
method [9] , we also evaluate it numerically by develop- 
ing an algorithm, the pair partitioning, inspired in the 
Trotter strategy used in quantum dynamics [10]. We 
split the kinetic terms to rewrite the Hamiltonian as 

H = (A.l) 

i 

Now, each term represents an effective Hamiltonian 
for two coupled oscillators, with twice the mass and 
half of the natural frequency each. Pair dynamics is 
solved analytically and impose a periodic evolution se- 
quence that alternates each coupled pair according to 
their parity. The total energy is not exactly conserved 
but fluctuates with an amplitude Ai5 around the ideal 
conserved value. Since AiJ is proportional to 5t^ , the 
square of the temporal step, it becomes negligible for 
typical cases where 5t = O.OIlj,^^- The fact that each 
piecelike dynamics is perfectly reversible is very impor- 
tant for the test of different time reversal procedures. 
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